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Excitation-contraction coupling (ECC) in tlie cardiac myocyte is mediated by a 
number of Iniglily integrated mechanisms of intracellular Ca^+ transport. Voltage- and 
Ca^+ -dependent L-type Ca^+ channels (LCCs) allow for Ca^+ entry into the myocyte, 
which then binds to nearby ryanodine receptors (RyRs) and triggers Ca^+ release from 
the sarcoplasmic reticulum in a process known as Ca^+-induced Ca^+ release. The highly 
coordinated Ca^+-mediated interaction between LCCs and RyRs is further regulated by 
the cardiac isoform of the Ca^+/calmodulin-dependent protein kinase (CaMKII). Because 
CaMKII targets and modulates the function of many ECC proteins, elucidation of its 
role in ECC and integrative cellular function is challenging and much insight has been 
gained through the use of detailed computational models. Multiscale models that can 
both reconstruct the detailed nature of local signaling events within the cardiac dyad 
and predict their functional consequences at the level of the whole cell have played an 
important role in advancing our understanding of CaMKII function in ECC. Here, we review 
experimentally based models of CaMKII function with a focus on LCC and RyR regulation, 
and the mechanistic insights that have been gained through their application. 
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INTRODUCTION 

Cardiac electrophysiology is a discipline with a rich and deep 
history dating back more than a half-century. Much of our 
understanding of fundamental biological mechanisms in this 
field comes from experimental research coupled with integrative 
mathematical modeling. The goal of this modeling has been to 
achieve a quantitative understanding of the functional mecha- 
nisms and relationships that span from the level of molecular 
structure and function to integrated cardiac myocyte behavior in 
health and disease. Advances in computational techniques and 
hardware have led to the recent development of many sophis- 
ticated large-scale models of heart tissue, in which these elec- 
trophysiological myocyte models are the fundamental building 
blocks. Some of the most fundamental advances in computa- 
tional cell biology, including formulation of dynamic models of 
voltage-gated ion channels, mitochondrial energy production, 
membrane transporters, intracellular Ca^+ dynamics, and sig- 
nal transduction pathways have emerged from this field. One 
such signaling pathway is that involving the Ca^+ /calmodulin 
(CaM)-dependent protein kinase (CaMKII), which targets and 
functionally regulates a number of proteins that are central to 
cardiac excitation-contraction coupling (ECC), and has been 
identified as a promising new target for antiarrhythmic therapy 
in patients with heart failure (Rokita and Anderson, 2012). This 
review focuses on the role of experimentally based models of 
CaMKII function with a focus on L-type Ca^+ channel (LCC) 
and ryanodine receptor (RyR) regulation, and the mechanistic 
insights that have been gained through their application. 



The nature of ECC is linked closely to the micro-anatomical 
structure of the cell. Sarcomeres, the basic units of contractile pro- 
teins, are bounded on both ends by the t-tubular system (Soeller 
and Cannell, 1999; Bers, 2001). The t-tubules extend deep into 
the cell and approach the sarcoplasmic reticulum (SR), an intra- 
cellular luminal organelle involved in the uptake, sequestration 
and release of Ca^+. The junctional SR (JSR) is the portion of 
the SR most closely approximating (within 12-15 nm Franzini- 
Armstrong et al., 1999) the t-tubules. The close proximity of 
these two structures forms a microdomain known as the dyad. 
RyRs are Ca^+ -sensitive Ca^+ -release channels which are pref- 
erentially located in the dyadic region of the JSR membrane. 
In addition, LCCs are preferentially located within the dyadic 
region of the t-tubules, where they are in close opposition to 
the RyRs. The process by which Ca^+ enters the myocyte dur- 
ing the initial depolarization stages of the action potential (AP) 
via voltage-gated LCCs and triggers the opening of RyRs, and 
hence JSR Ca^+ release, is known as Ca^+ -induced Ca^+ release 
(CICR) and comprises the initial transduction event of ECC. 
Recent imaging studies have demonstrated that dyadic clefts are 
small, where junctions contain an average of 14 RyRs, with some 
clusters incompletely filled with RyRs, and with a large fraction of 
clusters closely spaced within 20-50 nm, suggesting that smaller 
clusters may act together to function as a single site of CICR 
(Baddeley et al., 2009; Hayashi et al., 2009). 

The elucidation of CICR mechanisms became possible with 
the development of experimental techniques for simultaneous 
measurement of L-type Ca^+ current (IcaL) and Ca^+ transients, 
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and detection of local Ca^+ release events known as Ca^+ sparks 
(Cheng et al., 1993). The evidence for tight regulation of SR Ca^^ 
release by triggering LCC current (for a review see Soeller and 
Cannell, 2004) gave rise to and later verified the local control 
theory of ECC (Stern, 1992). The mechanism of local control 
predicts that tight regulation of CICR is achieved because LCCs 
and RyRs are regulated by local dyadic Ca^+. It is now widely 
accepted that CaMKII, a ubiquitous Ca^+ -dependent protein that 
can become highly activated in the dyad (Currie et al., 2004; 
Hudmon et al., 2005; Maier, 2005), can modulate CICR via phos- 
phorylation of a number of ECC proteins including LCCs and 
RyRs (Maier and Bers, 2007). CaMKII modulates LCCs via a 
Ca^+ -dependent positive-feedback regulatory mechanism known 
as IcaL facilitation (Anderson et al, 1994; Xiao et al, 1994; 
Yuan and Bers, 1994). This is observed as a positive staircase 
in current amplitude in combination with a slower rate of IcaL 
inactivation upon repeated membrane depolarization. Although 
multiple mechanisms have been suggested to underlie IcaL facili- 
tation, experiments from Dzhura et al. (2000) demonstrated that 
CaMKII-dependent IcaL facilitation alters the behavior of LCCs 
such that high activity gating modes with prolonged open times 
are more likely to occur. The specific functional effects of CaMKII 
phosphorylation of RyRs remain controversial. CaMKII has been 
shown to increase (e.g., Wehrens et al., 2004) or decrease (e.g., 
Lokuta et al., 1995) the channel open probability. Inconsistencies 
in experimental findings regarding the role of CaMKII on RyR 
function may not be surprising given the methodological differ- 
ences between the studies. In addition, the high sensitivity of RyR 
gating to both cytosolic and SR Ca^+ makes it difficult to interpret 
experiments in which these Ca^+ concentrations are not tightly 
controlled. In one well designed study, Guo et al. (2006) devel- 
oped an experimental protocol that overcame this challenge and 
provided strong evidence that CaMKII mediated phosphoryla- 
tion of RyRs increases channel open probability. Such an increase 
is also consistent with transgenic studies of CaMKII overexpres- 
sion which report increased Ca^^ spark frequency in response 
to elevated CaMKII activity (Maier et al., 2003; Kohlhaas et al, 
2006). 

The CaMKII holoenzyme exists as an elaborate macromolec- 
ular complex consisting of two stacked ring-shaped hexamers 
(Hudmon and Schulman, 2002). Each of the holoenzyme's 12 
subunits can be activated through the binding of Ca^+ -bound 
CaM (Ca^+/CaM) to the CaMKII regulatory domain in response 
to beat-to-beat transient increases of intracellular Ca^+ con- 
centration. In cardiac myocytes, an activated CaMKII molecule 
can undergo autophosphorylation by neighboring subunits at 
threonine amino acid residues in its regulatory domain, which 
allows the kinase to retain activity even upon dissociation of 
Ca^+/CaM (Meyer et al., 1992). An alternative mechanism of 
oxidative CaMKII activation has also been identified (Erickson 
et al, 2008). As a result of the fact that CaMKII regulates multi- 
ple protein targets (both directly and indirectly involved in ECC), 
and that the effect of CaMKII mediated phosphorylation on any 
particular target protein may involve complex alterations in bio- 
physical function (e.g., mode-switching behavior of LCCs), the 
task of elucidating its role in cardiac function in both normal and 
failing hearts continues to present great challenges and cannot 
be accomplished via experiments alone. Some recent advances in 



understanding the mechanisms of CaMKII-dependent function 
in cardiac myocytes have been discovered by coupling mathemat- 
ical models to experimental observations. In this review we will 
focus on key models of CaMKII-mediated regulation of LCCs 
and RyRs arising from both Ca^+/CaM-dependent as well as 
oxidative activation pathways, the integration of such models into 
whole-cell models, and the mechanistic insights that have been 
obtained from this work. A summary of the models covered here 
is provided in Table 1 . 

MODELS OF Ca^+/CaM-DEPENDENT CaMKII ACTIVATION 
AND REGULATION 

The first model of CaMKII signaling within the context of ECC 
integrated into the cardiac myocyte was presented by Hund and 
Rudy (2004). This model, known as the HRd model, incorpo- 
rates a scheme based on the work of Hanson et al. (1994), where 
a single population of CaMKII transitions from an inactive to 
active state in response to elevated subspace Ca^+ levels. CaMKII 
activity is assumed to modify the function of LCCs and RyRs by 
slowing their inactivation kinetics, thereby enhancing IcaL and 
JSR Ca^+ release, respectively. SR Ca^+ uptake was also enhanced 
by CaMKII activity via its action on phospholamban (PLB) and 
the SR Ca2+-ATPase (SERCA). The model predicted that CaMKII 
plays a critical role in the rate-dependent increase of the cytoso- 
lic Ca^+ transient by increasing ECC gain, but that it does not 
play a significant role in rate-dependent changes of AP duration. 
Similar results were obtained in the human cardiac AP model 
of O'hara et al. (2011), which incorporates the role of CaMKII 
phosphorylation on LCCs, the fast Na+ current (iNa)i and the 
transient outward K+ current (Itoi). The HRd model was later 
updated and used to investigate the role of altered CaMKII sig- 
naling in the canine infarct border zone (Hund et al., 2008). 
Experimentally measured elevation of CaMKII autophosphory- 
lation in this region is reproduced by the model, which indicates 
that hyperactive CaMKII impairs Ca^+ homeostasis by increasing 
Ca^+ leak from the SR. A number of additional modeling studies 
have focused on the role of CaMKII on various aspects of intracel- 
lular Ca^+ dynamics. Iribe et al. (2006) incorporated a conceptual 
CaMKII model into their previous myocyte model (Noble et al., 
1991) combined with a model of mechanics (Rice et al., 1999) in 
order to investigate its role on SR Ca^+ handling and interval- 
force relations. They found that a relatively slow time-dependent 
inactivation of CaMKII allowed for the reconstruction of a vari- 
ety of interval-force relations, including alternans. Koivumaki 
et al. (2009) built a model of the murine cardiac myocyte to 
analyze genetically engineered heart models in which CaMKII- 
mediated phosphorylation of LCCs is disrupted or CaMKII is 
overexpressed. They demonstrated how these genetic manipula- 
tions lead to the observed experimental phenotypes as a result 
of autoregulatory mechanisms that are inherent in intracellular 
Ca^+ cycling (e.g., steady-state regulation of SR content via Ca^+ 
release dependent inactivation of LCCs), and that disruption of 
the regulatory system itself (e.g., via CaMKII overexpression) 
leads to the most aberrant physiological phenotypes. Livshitz and 
Rudy (2007) formulated a new model of SR Ca^+ release kinet- 
ics and incorporated it into the HRd model in order to better 
understand regulation of Ca^+ and electrical alternans under var- 
ious pacing protocols. They found that increased CaMKII activity 
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Table 1 | Key features of CaMKII models in cardiac ECC. 



Cell model 

Hund and Rudy, 2004 
Iribe etal., 2006 
Livshitz and Rudy, 2007 

Grandi et al., 2007 

Saucerman and Bers, 2008 

Hund etal., 2008 
Koivumaki etal., 2009 
Hashambhoyetal., 2009 

Christensen etal., 2009 
Soltis and Saucerman, 2010 

Hashambhoyetal., 2010 
O'hara et al., 2011 

Hashambhoy et al., 2011 
Zang etal., 2013 
Morotti etal., 2014 



Species 

Dog 

Guinea pig 
Guinea pig/Dog 

Rabbit 

Rabbit 

Dog 

Mouse 
Dog 

Dog 
Rabbit 

Dog 
Human 

Dog 
Dog 

Mouse 



Parent ECC model 

Lug and Rudy, 1994 
Noble etal., 1991 
Luo and Rudy, 1994; Hund 
and Rudy, 2004 
Shannon et a!., 2004 

Shannon et a!., 2004 

Hund and Rudy, 2004 
Bondarenkoetal., 2004 
Greenstein and Winslow, 
2002 

Hund et al., 2008 
Saucerman and Bers, 2008 

Hashambhoyetal., 2009 
Hund and Rudy, 2004; Decker 
etal., 2009 

Hashambhoy et al., 2010 
Hund and Rudy, 2004 
Soltis and Saucerman, 2010 



CaMKII activation model 

Hanson et al., 1994 
Conceptual kinetic model 
Hanson et al., 1994 

Static formulation 
Target parameters adjusted 
Adapted from Dupont et al. 
(2003), Gaertner et al. (2004) 
Hanson et al., 1994 
Bhalla and Iyengar, 1999 
Dupont etal., 2003 

Hanson et al., 1994 

Dupont et a!., 2003; Gaertner 

etal., 2004 

Dupont et a!., 2003 

Hund and Rudy, 2004 

Dupont etal., 2003 
Hanson et al., 1994 
Dupont et al., 2003; Gaertner 
etal., 2004 



Included CaMKII targets 

LCCs, RyRs, SERCA/PLB 
LCCs, RyRs, SERCA 
LCCs, RyRs, SERCA/PLB 

ICaL. iNa, Ito 

no target 

LCCs, RyRs, SERCA/PLB, iNa 
LCCs, RyRs, SERCA/PLB 
LCCs 

LCCs, RyRs, SERCA/PLB, iNa 
LCCs, RyRs, PLB, iNa, Ito 

LCCs, RyRs, PLB 

LCCs, PLB, iNa, Ito 

LCCs, RyRs, PLB, Ins 

LCCs, RyRs, iNa, Ito, ki. SERCA 

LCCs, RyRs, PLB, Ins, Ito-ki, Incx 



leads to increased alternans magnitude as well as the appearance 
of both Ca^^" and electrical alternans at lower pacing rates (where 
they would not normally occur) . This model identifies Ca^+ alter- 
nans as the underlying mechanism for electrical alternans, both of 
which were eliminated with CaMKII inhibition, suggesting this as 
an antiarrhythmic strategy. Zang et al. (2013) recently developed 
a new canine myocyte model, also based on the HRd model, to 
study the role of upregulated CaMKII in heart failure. Similarly, 
they find that enhanced RyR Ca^+ -sensitivity and SR Ca^+ leak 
mediated by CaMKII overexpression alters Ca^+ handling in a 
manner that promotes alternans, while AP prolongation occurs 
primarily due to an associated down regulation of K+ currents. 
Interestingly, they find that blocking SR Ca^+ leak restores con- 
traction and relaxation function, but does not eliminate alternans 
completely. 

Grandi et al. (2007) developed a model of CaMKII overexpres- 
sion in the rabbit ventricular myocyte, which incorporates the 
functional effects of CaMKII-mediated phosphorylation of LCCs, 
as well as I^a and Itoi- This model shows that while CaMKII- 
mediated action on LCCs prolongs the AP, the combined effect of 
CaMKII on all three of the targets studied leads to AP shortening. 
While this study was primarily motivated by and focused on 
the functional role of CaMKII-mediated regulation of iNa (see 
accompanying article in this series by Grandi and Herren), 
it clearly demonstrates that the multiple targets of CaMKII 
interact to produce net alterations of integrated cellular function 
which are difficult to predict from knowledge if its functional 
effects on isolated individual targets. Saucerman and Bers (2008) 
developed and incorporated models of CaM, CaMKII, and 
calcineurin (CaN) into the rabbit ventricular myocyte model of 



Shannon et al. (2004) in order to better understand the functional 
consequences of the different affinities of CaM for CaMKII and 
CaN during APs. The model predicts that the relatively high Ca^+ 
levels that are achieved in the cardiac dyad lead to a high degree 
of CaM activity which results in frequency-dependent CaMKII 
activation and constitutive CaN activation, whereas the lower 
Ca^+ levels in the cytosol only minimally activate CaM, which 
allows for gradual CaN activation, but no significant activation of 
CaMKII. The prediction that robust beat-to-beat oscillations of 
CaMKII activity occur in the dyad (i.e. in the vicinity of RyRs and 
LCCs) but not in the cytosol is a key factor that would influence 
the way in which local signaling mechanisms would be incor- 
porated into models that followed that included detailed recon- 
struction of dyadic Ca^^" dynamics. In one such model, Soltis 
and Saucerman (2010) integrated dynamic CaMKII-dependent 
regulation of LCCs, RyRs, and PLB with models of cardiac ECC, 
CaMKII activation, and /j-adrenergic activation of protein kinase 
A (PKA). In this model phosphorylation of all CaMKII substrates 
exhibits positive frequency dependence similar to that observed 
in experiments (De Koninck and Schulman, 1998). However, 
both CaMKII activity and target protein phosphorylation levels 
adapt to changes in pacing rate with distinct kinetics (see Figure 
2 of Soltis and Saucerman, 2010). CaMKII-mediated phosphory- 
lation is relatively fast at LCCs, slower at RyRs, and very slow at 
PLB. Additionally, the model predicts a high degree of phospho- 
rylation of LCCs in control conditions, but a moderate amount 
of phosphorylation of RyRs and PLB (e.g., PLB phosphorylation 
<10% at all frequencies). In addition, this study predicts a novel 
mechanism in which CaMKII and PKA synergize to form a 
positive feedback loop of CaMKII-Ca2+ -CaMKII regeneration. 
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Furthermore, this model predicts that CaMKII-mediated hyper- 
phosphorylation of RyRs, which renders them leaky, may be 
a proarrhythmic trigger via induction of delayed afterdepolar- 
izations (DADs). Recently, Morotti et al. (2014) modified the 
Soltis and Saucerman (2010) rabbit model to farther investigate 
the arrhythmogenic role of another synergistic interaction 
in mouse ventricular myocytes. This synergism is that of the 
positive feedback loop of CaMKII-Na+-Ca^+ -CaMKII in which 
CaMKII-dependent increases in intracellular Na+ level perturbs 
Ca^+ homeostasis and CaMKII activation. Simulation results 
from Morotti et al. (2014) demonstrate that the feedback between 
disrupted Na^" fluxes and CaMKII signaling is exaggerated when 
CaMKII is overexpressed. Under this condition, and upon an 
increase in intracellular Na+ concentration,the model predicts 
Ca^+ overload and enhancement of CaMKII activity which in 
turn increases spontaneous Ca^+ release events via an increase in 
RyR phosphorylation. This CaMKII-dependent hyperphospho- 
rylation of RyRs exacerbates the associated electrophysiological 
instability as simulated by the occurrence of DADs (see Figure 
7 of Morotti et al, 2014). In this model, the DADs do not 
occur when CaMKII target phosphorylation is clamped to basal 
levels. 

Recently, Hashambhoy et al. (2009) developed a biochemically 
and biophysicaUy detailed model of CaMKII-mediated phospho- 
rylation of LCCs. The model includes descriptions of the dynamic 
interactions among CaMKII, LCCs, and phosphatases as a func- 
tion of dyadic Ca^+ and CaM levels, and has been incorporated 
into an integrative model of the canine ventricular myocyte with 
stochastic simulation of LCC and RyR channel gating within a 
population of release sites based on the theory of local control of 
ECC (Greenstein and Winslow, 2002). A schematic representa- 
tion of one such model Ca^+ -release site is shown in Figure lA. 
This CaMKII-LCC model is formed by the integration of three 
modules: a CaMKII activity model which reflects the different 
structural and functional states of the kinase based on the work 
of Hudmon and Schulman (2002) and Dupont et al. (2003), an 
LCC phosphorylation model derived from studies of facilitation 
in LCC mutants (Grueter et al., 2006; Lee et al, 2006), and a pre- 
viously developed LCC gating model (Jafri et al, 1998; Greenstein 
and Winslow, 2002). In this model it is assumed that there is one 
12-subunit CaMKII holoenzyme tethered to each LCC (Hudmon 
et al, 2005), each CaMKII monomer can transition among a 
variety of activity states (see Figure 2 of Hashambhoy et al., 
2009), and CaMKII monomers can catalyze phosphorylation of 
individual LCCs. The LCC gating model reflects two forms of 
channel gating, mode 1 (normal activity) and mode 2 (high activ- 
ity with long openings). In the model, LCC phosphorylation 
promotes transitions of LCCs from mode 1 to mode 2 gating, 
based on experimental observations with constitutively active 
CaMKII (Dzhura et al, 2000). This model demonstrates that 
these CaMKII-dependent shifts of LCC gating patterns into high 
activity gating modes may be the underlying mechanism of IcaL 
facilitation. In addition, the model predicts that this CaMKII- 
mediated shift in LCC gating leads to an apparent increase in 
the speed of both IcaL inactivation and recovery from inacti- 
vation, both of which are experimentally associated with IcaL 
facilitation, with no change to the underlying intrinsic LCC 



inactivation kinetics. Hashambhoy et al. (2010) further expanded 
this model to include CaMKII-dependent regulation of RyRs. In 
this model update, RyR phosphorylation is modeled as a function 
of dyadic CaMKII activity and it is assumed that RyR sensitiv- 
ity to dyadic Ca^+ is increased upon phosphorylation, the degree 
to which is constrained by experimental measurements of Ca^+ 
spark frequency and steady state RyR phosphorylation (Guo et al., 
2006). This study demonstrated that under physiological condi- 
tions, CaMKII phosphorylation of LCCs ultimately has a greater 
effect on ECC gain, RyR leak flux, and AP duration as com- 
pared with phosphorylation of RyRs (Figure IB). AP duration 
at 90% repolarization (APD90) correlates well with a CaMKII- 
mediated shift in modal gating of LCCs (Figure IC). A modest 
additional increase in LCC phosphorylation, and hence mode 2 
gating, beyond that shown in Figure IC leads to the appearance 
of early afterdepolarizations (EADs) in simulated APs (see Figure 
5 of Hashambhoy et al., 2010). The results of these model anal- 
yses suggest that LCC phosphorylation sites may in fact prove to 
be a more effective target than those on the RyR for modulat- 
ing diastolic RyR-mediated Ca^+ leak and preventing abnormal 
proarrhythmic cellular depolarization. 

MODELS OF OXIDATIVE CaMKII ACTIVATION AND 
REGULATION 

The modeling studies described above focus on the role 
of CaMKII-mediated regulation of its targets via the classi- 
cal Ca^+/CaM-dependent activation pathway which involves 
CaMKII autophosphorylation. Recently a novel mechanism 
for oxidative CaMKII activation was discovered in the heart 
(Erickson et al., 2008). This alternative oxidation-dependent 
pathway involves the oxidation of CaMKII at specific methio- 
nine residues, which has been shown to produce persistent 
kinase activity and increase the likelihood of arrhythmias (Xie 
et al, 2009; Wagner et al., 2011). These new findings impli- 
cate oxidative CaMKII activation as a putative mechanistic link 
between the accumulation of reactive oxygen species (ROS) 
and life-threatening cardiac arrhythmias (Erickson et al., 2011). 
Therefore, developing a quantitative understanding of the role 
of CaMKII oxidation in regulating cardiac ECC requires the 
integration of computational models that link cellular ROS and 
redox balance, CaMKII activity and function, ECC, and whole- 
cell electrophysiology. Such models will prove to be powerful tools 
for teasing out important mechanisms of arrhythmia in heart 
disease. 

In one recent study, Christensen et al. (2009) developed a 
deterministic model of oxidative CaMKII activation based on 
their own experimental findings to study the role of this signaling 
pathway in the canine infarct border zone following myocar- 
dial infarction. A new simplified model of CaMKII activation 
was developed based on the work of Dupont et al. (2003), and 
incorporated into a model of a single cardiac fiber. Each myocyte 
within the fiber model was represented by the Hund et al. (2008) 
model, which as noted above, incorporates CaMKII effects on lNa> 
IcaL. SR Ca^+ leak and uptake, as well as ion channel remod- 
eling associated with the infarct border zone. Simulation results 
demonstrate that enhanced oxidative CaMKII activity is asso- 
ciated with reduced conduction velocity, increased refractory 
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FIGURE 1 I (A) Schematic representation of the cardiac dyad. Ca^+ ions 
pass through the LCC and enter the dyad, where they bind to RyRs, 
triggering CICR. The activity of each CaMKII monomer is a function of 
dyadic Ca^+ levels. Ca^+ ions released from the SR diffuse out of the dyad 
into the bulk myoplasm. (B) Simulated results from a 1-Hz AP-pacing 
protocol under different LCC and/or RyR phosphorylation conditions. 
(C) APDgo as a function of average LCC phosphorylation levels. Fully 
phosphorylated LCCs gate in mode 2, which exibit long-duration openings. 
(A) Reprinted with permission from Hashambhoy et al. (2009). Copyright 
2009 Elsevier. (B,C) Reprinted with permission from Hashambhoy et al. 
(2010). Copyright 2010 Elsevier. 



periods, and a greater likelihood of conduction block formation. 
Furthermore, CaMKII inhibition in the model improves conduc- 
tion and reduces vulnerability to conduction block and reentry. 
The results of Christensen et al. (2009) are attributed primarily 
to CaMKII-mediated regulation of iNa kinetics and availability. 
They note that CaMKII activation will also impact ECC proteins 
in ways that may promote arrhythmias, but these mechanisms 
were not explored in detail in this model. 

Along these lines, Foteinou et al. (2013) have recently devel- 
oped a novel stochastic model of CaMKII activation (Figure 2A) 
that reflects the functional properties of the cardiac isoform 
including both the canonical phosphorylation-dependent acti- 
vation pathway and the newly identified oxidation-dependent 
pathway. This model builds upon the work of Hashambhoy et al. 
(2010, 2011) in order to incorporate recent experimental data for 
CaM affinity and autophosphorylation/oxidation rates measured 
specifically for the cardiac isoform of CaMKII (Gaertner et al, 
2004; Erickson et al., 2008). Predicated upon this, the authors 
implemented the four-state deterministic activation model of 
Chiba et al. (2008) within a stochastic framework that is con- 
strained by the geometry of the CaMKII holoenzyme. This was 
accomplished by restricting CaMKII autophoshorylation events 
to occur only between adjacent CaMKII subunits. The stochas- 
tic activation model was further modified by including oxidized 
active states in addition to a Ca^^/CaM-bound active state, an 
autophosphorylated Ca^+/CaM-bound state, and an autophos- 
phorylated Ca^+/CaM-dissociated state (i.e., an autonomous 
active state). Following incorporation into the canine cardiac 
myocyte model (Hashambhoy et al., 2011), Foteinou et al. (2013) 
obtain steady state APs at a pacing cycle length (PCL) of 2 s 
as illustrated in Figure 2B. This model predicts increased like- 
Hhood of EADs under this protocol upon increased oxidative 
stress (application of 200 |iM H2O2, which is at the high end of 
pathophysiological levels), as shown in Figure 2C. Notably, the 
model simulates no EADs in the presence of ROS at a PCL of 
1 s, which is consistent with the experimentally measured rate- 
dependence of EADs in the presence of 200 [xM H2O2 (Zhao 
et al., 2012). The model predicts that EADs result from both 
enhanced CaMKII-dependent shifts of LCCs into highly active 
mode 2 gating and a progressively increased late component 
of Na+ current (iNah)- The simulated increase in IcaL via oxi- 
dized CaMKII corroborates the experimental findings of Song 
et al. (2010) who demonstrated that oxidative CaMKII activa- 
tion is involved in the regulation of LCCs. In particular, the in 
silico model of H2O2 exposure discussed here predicts that the 
fraction of LCCs gaiting in mode 2 will shift from 5% in con- 
trol (absence of ROS) to 7% upon this increase in ROS. Using 
this model, the simulated diastolic RyR leak also increases in 
response to oxidative stress (~44%) primarily as an indirect con- 
sequence of increased Ca^+ influx via LCCs. Wagner et al. (201 1) 
recently observed a ~15-fold H202-mediated increase in SR 
Ca^+ leak, which is far greater than that produced in this model 
(<2-fold). They, however, provide additional evidence indicating 
that their observed increase in Ca^+ leak does not require the 
presence of CaMKII, suggesting an important role for CaMKII- 
independent mechanisms of ROS-mediated alteration of cardiac 
ECC as well. 
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FIGURE 2 I (A) State diagram of the stochastic CaMKII activation model of 
Foteinou et al. (2013). Prior to the introduction of Ca^+/CaM, all of the CaMKII 
subunits are in the inactive form (state I). Activation occurs upon binding of 
Ca^+/CaM (state B), autophosphorylation (state P), and oxidation (state Oxb). 
Autonomous active states (Ca^+/CaM-unbound) can be either 
autophosphorylated (state A) or oxidized (state Oxa). The model also includes 



an active state that is both oxidized and phosphorylated (state Oxp). 

(B) Simulated steady state APs under control conditions (0|iM H2O2). 

(C) Simulated APs, some of which exhibit EADs, under conditions of elevated 
oxidant stress (200 |j,M H2O2). Simulated APs are from a 2 s PCL pacing 
protocol (ensemble of 12500 calcium release units). For each condition, 
results for 10 consecutive APs following an initial 10s of pacing are shown. 



As these ROS-mediated effects need to be farther exam- 
ined, future studies will focus on establishing quantitative 
links between cellular ROS and redox balance, CaMKII activity 
and function, ECC, and whole-cell electrophysiology. Recently, 
Gauthier et al. (2013a,b) developed mechanistic models of ROS 
production and scavenging to investigate how these two com- 
peting processes control ROS levels in cardiac mitochondria. 
Simulations confirm the hypothesis that mitochondrial Ca^+ 
mismanagement leads to decreased scavenging resources and 
accounts for ROS overflow as is believed to occur in heart failure 
(Hill and Singal, 1996). Notably the ROS regulation module of 
Gauthier et al. (2013b) enables its use in larger scale heart mod- 
els designed to simulate and study how mitochondrial ROS and 
the functional consequences of its accumulation, such as CaMKII 
oxidation, regulate cellular physiological function, AP properties, 
and arrhythmogenesis. 

CONCLUSION 

Integrative modeling of cardiac ECC, cell signaling, and myocyte 
physiology has played a critical role in revealing mechanis- 
tic insights across a range of biological scales. With respect to 
CaMKII function at the smallest scale, models have shed light 
on how Ca^+ ions, CaM, CaMKII, LCCs, and RyRs interact in 
the cardiac dyadic junction both at rest and during triggered ECC 
events. On an intermediate scale, models have predicted the con- 
sequences of normal and abnormal CaMKII signaling on whole- 
cell Ca^+ cycling, effects of ROS imbalance, AP shape, and the 
generation of cellular arrhythmias such as EADs. Incorporation of 



these cellular models into higher scale tissue simulations has pro- 
vided important insight into the relationship between CaMKII 
function, electrical wave conduction velocity, and the emergence 
of arrhythmogenic substrates in diseased tissue. The continuum 
of biological scales spanned by these models allows for the devel- 
opment of multiscale approaches whereby we can predict and 
understand the emergence of macroscale phenotypes as a conse- 
quence of CaMKII-mediated molecular signaling events. A great 
deal of evidence now implicates CaMKII as a nexus point linking 
heart failure and arrhythmias (Swaminathan et al, 2012), iden- 
tifying it as a prime target for antiarrhythmic therapies (Fischer 
et al, 2013). Despite this growing body of evidence, the model- 
ing studies presented here demonstrate that it remains difficult 
to identify which of the many CaMKII target proteins are pri- 
marily responsible for the functional changes that increase the 
likelihood of arrhythmia development. Some studies (Hund et al., 
2008; Soltis and Saucerman, 2010; Zang et al, 2013; Morotti 
et al., 2014) identify phosphorylation of RyRs, which leads to ele- 
vated SR Ca^+ leak and/or increased likelihood of spontaneous 
Ca^+ release events and DADs, as the key mechanism underlying 
arrhythmogenesis. Others (Hashambhoy et al, 2010; Foteinou 
et al., 2013) predict that LCC phosphorylation, which leads to ele- 
vated channel open probability via high activity gating, plays the 
predominant role by altering whole cell Ca^+ levels and AP prop- 
erties, including the appearance of EADs. These mechanisms may 
co-exist in heart failure, and it is likely that both play an impor- 
tant role. Moreover, CaMKII modulation of other targets such as 
PLB, iNa, and K+ currents, many of which have been incorporated 
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into the above cell models, adds additional layers of complexity to 
the task of interpreting and predicting the mechanisms by which 
CaMKII signaling alters cardiac function in health and disease. 
As a result of the breadth and diversity of CaMKII targets, strate- 
gies that involve broad inhibition of CaMKII activity will impact 
and alter the function of many cellular subsystems, with com- 
plex consequences that may not all be beneficial. As mathematical 
models of CaMKII signaling in the cardiac dyad and beyond are 
further developed, they will play a key role in identifying target- 
specific therapeutic strategies for improving myocardial function 
and reducing arrhythmias. 
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